Anderson‐accelerated polarization schemes for fast Fourier transform‐based computational homogenization

Classical solution methods in fast Fourier transform-based computational micromechanics operate on, either, compatible strain fields or equilibrated stress fields. By contrast, polarization schemes are primal-dual whose iterates neither nor equilibrated. Recently, it was demonstrated that may outperform the classical methods. Unfortunately, their power critically depends on a judicious choice of numerical parameters. In this work, we investigate extension by Anderson acceleration and demonstrate combination leads to robust general-purpose solvers for micromechanics. We discuss (theoretically) optimum parameter methods, describe how fits into picture, exhibit characteristics newly designed problems industrial scale interest.